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Abstract 

It is well known that the building blocks for state sum models of quantum gravity are given 
by 6j and lOj symbols. In this work we study the asymptotics of these symbols by using their 
expressions as group integrals. We carefully describe the measure involved in terms of invariant 
variables and develop new technics in order to study their asymptotics. Using these technics we 
compute the asymptotics of the various Euclidean and Lorentzian 6j-symbols. Finally we compute 
the asymptotic expansion of the lOj symbol which is shown to be non-oscillating, in agreement 
with a recent result of Baez et al. We discuss the physical origin of this behavior and a way to 
modify the Barrett-Crane model in order to cure this disease. 
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I. INTRODUCTION 



It is well known that state sum models defined over a triangulation of space time are a 
central tool in the construction of discrete transition amplitudes for quantum gravity. It was 
first understood by Ponzano and Regge [1] that the transition amplitude of 3d Euclidean 
gravity could be expressed as a state sum model defined over a triangulation of space time. 
This state sum model is obtained by summing, over SU(2) representations labelling the 
edges of the triangulation, a given weight depending on this labelling. This sum over repre- 
sentations is interpreted as a sum over geometry, and the weight is a product of SU(2) 6j 
symbols associated with each colored tetrahedron. This construction can be extended to the 
case of 4d gravity and 2 + 1 Lorentzian gravity. It was shown in [2, 3] that the building block 
in the construction of transition amplitudes for 2 + 1 Lorentzian gravity are the SL(2, M) 6j 
symbols. In the context of 4 dimensional gravity, Barrett and Crane [4] proposed to use the 
so-called lOj symbol in order to construct the discretized transition amplitudes. 

One of the key argument in favor of the Ponzano- Regge model as a model for 3d gravity 
is the fact that the asymptotic behavior of the 6j symbol reproduces the discretized Regge 
action for 3d gravity. This asymptotic has been conjectured in 1968 by Ponzano and Regge 
but proved only in 1999 by Roberts [5]. In the same way, one of the argument in favor of 
the Barrett-Crane amplitude as a building block for quantum gravity was the fact shown 
by Barrett and Williams [6] that the stationary contribution to the asymptotic of the lOj 
symbol reproduces also the discretized Regge action of 4d gravity. However, recent numerical 
simulations and computations by Baez et al. [7, 8] have shown that the asymptotic of the 
lOj symbol is not dominated by an oscillating contribution. 

In this paper we address this issues by using expressions for the 6j and lOj symbol as 
integrals over group elements [9]. This is done by carefully describing the measure involved 
in terms of gauge invariant variables. We show that the asymptotics separates into a con- 
tribution coming from a stationary phase approximation and a contribution coming from 
degenerates configurations associated with singularities of the integrand. The stationary 
phase produces an oscillatory Regge behavior, whereas the degenerate contribution is non 
oscillating. It should be noted that this general method has been outlined in the recent 
paper of Baez, Christensen and Egan [8], and even if our approach is independent from 
their, one can view our paper as giving proofs of the conjectures they made. Also, while we 
were completing the redaction of this work, we become aware of the very recent article [10], 
which is leading to the same conclusions as ours but using a different approach. 

In section II we express the square of the SU(2) 6j-symbol as an integral over the space of 
spherical tetrahedra. We show that the measure of integration is given by the inverse square 
root of the determinant of the Graham matrix [11] associated with spherical tetrahedra. We 
show that the integral naturally separates into two parts : one for which the asymptotics 
is dominated by an oscillating contribution, associated with fiat non-degenerate tetrahedra 
and obtained by a stationary phase approximation; and one for which the asymptotics 
comes from boundary contributions which label degenerate tetrahedra. We prove that this 
contribution is given by an integral associated with the Euclidean group, as was conjectured 
by Baez et al [8]. 

In section III we express the square of the SL(2,R) 6j-symbol as an integral over the 
space of AdS tetrahedra. We also split this integral in two parts, and show that the station- 
ary phase is dominated by an oscillating contribution associated with flat non-degenerate 
Lorentzian tetrahedra. The phase of this contribution is the Lorentzian Regge action and the 
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module is the volume of the Lorentzian tetrahedron. In section IV we apply our technics to 
the lOj-symbol and show that the leading contribution is a non-oscillating one, dominating 
the oscillating Regge action term found by Barrett and Williams. Finally in section V we 
discuss the physical meaning and origin of this non-oscillating behavior for the lOj-symbol, 
in the spirit of the statistical mechanics models of 'order by disorder', and we propose a 
modification of the Barrett-Crane model to avoid this problem and recover an oscillatory 
Regge action behavior. 



II. ASYMPTOTICS OF THE 6J-SYMBOL 

A. Integral expression for the square of the 6j-symbol 

We are interested in the computation of the asymptotics of the square of the 6j-symbol. 
The 6j-symbol is a real number which is associated to the labelling of the edges of a tetra- 
hedra by SU (2) representations. It is obtain by combining four normalized Clebsch-Gordan 
coefficients along the six edges of a tetrahedra (see [5]). Let us denote by Vi the SU{2) 
representation of spin 1/2 and by Xiio) = try; ((7) the character of this representation. Lets 
/ = 0, 1,2,3 be the 4 vertices of a tetrahedra T and (/J) the edges of T. We associate a 
representation Vijj to each edge (/J) of T. It is well known [9, 16] that the square of the 
6j-symbol can be expressed as the following Feynman integral over SU(2) 

likj) = I Y= [ [dgr] n Xi,A9J97'). (1) 

hs ks h2 j Jqa 

where the normalized Haar measure dg is used. Such an identity is clear since the integral 
over each group clement produces a pair of Clebsch-Gordan coefficient which are then com- 
bined into one tetrahedron and its mirror image. It is also clear from this expression that 
the /(//j) 7^ only if the \Iik — < ^/j < Iik + Ul and Ijk + Iil + hj is an even integer. 
In this expression, the integral is over four copies of the group, however the integrand is 
invariant under the transformations 

gi ^ hkgjk-\ (2) 

where h,k are SU{2) elements. It is therefore possible to gauge out this symmetry and 
write the integral purely in term of gauge invariant variables. A natural choice for the gauge 
invariant variables is given by the 6 angles 6ij e [0, tt]: 

cos Oij = ]^trv^ {gjgj^) ■ (3) 

Note that this gauge symmetry has a geometrical interpretation. The four group elements 
gi define 4 points, hence a tetrahedron, in S^. The angles 9ij are the spherical lengths of 
the edges of this tetrahedron, which indeed parameterize its invariant geometry. Notice that 
the symmetry includes refiexions g —g. The integral (1) can be written in terms of these 
angles as follows: 



Theorem 1 



/(M=|/in''M "-;:'f":''^"' (4) 

Jv^-^^ A/det cos^/j 
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where V^^ is the subset of [0; tt]^ of angles satisfying the relations : 

Oij < Oik + ejK, (5) 
27r > Oij + eiK + ejK, (6) 

for any triple (/, J, K) of distinct elements. Geometrically this domain is the set of all 

possible spherical tetrahedra. 

The denominator of this integral is the square root of the determinant of the 4x4 Graham 
matrix Gu = cos 9ij associated with the corresponding spherical tetrahedron . This deter- 
minant is zero if and only if 9ij belongs to the boundary of the domain D^, which is the set 
of degenerate tetrahedra having zero volume. Note that if we consider Ojj to be the dihedral 
angles of a spherical tetrahedra, dV^^ is the set of flat tetrahedra (see the remark 3 below) . 

Proof : — 

Notations and spherical geometry: The invariant measure is obtained by a Faddev-Popov 
procedure. In order to present the details of this procedure, we first need to do a httle bit of 

spherical geometry and introduce a new set of angles denoted 6i, Oij, 9ij e [0, tt], i = 1, 2, 3. 
This angles are related in the following way to the angle 9ij of eq.(3). 

— ^Oi (7) 
cos Oij — cos 9i cos 9j + sin 9i sin 9j cos 9ij , (8) 

cos 9ij — cos 9ik cos 9kj + sin Oik sin 9kj cos 9ij, (9) 

where in the eq.(9), i,j,k is any permutation of 1,2,3. The geometrical interpretation of 
these angles is as follows : first 9i are the lengths of the edges Oi, 9ij is the angle between 

the edges (Oi) and (Oj), finally 9ij is the dihedral angle of the edge (Ok) opposite to the edge 
{ij) (see fig.l). In order to see this, let (Tj, i = 1, 2, 3 be the Pauli matrices, aiaj = Sij+ie^jUk 
and let us introduce the normalized vectors fij e S"^ such that 

9i9o^ = e'^*"^-" = cos 9i + i sin 9i{ni ■ a). (10) 

Then 

cos 9ij = ^tr{gjg:r^) = ^tr {gjgo^{gigo^)~^) = cos 9i cos 9 j + sin 9isin 9 j{ni ■ fij), (11) 

so cos 9ij = Hi ■ fij. Let us introduce the unit normal vector to the face (Oij): = Ui A 
nj/\ni A fij I, where ijk is a cyclic permutation of 123. A simple computation shows that 

cos9ij — —tti ■ Qj, (12) 

which characterizes 9ij as the dihedral angle of the edge Ok. 

Measure: Using these geometrical elements, we can compute the invariant measure. In order 
to get this measure, we first have to fix the symmetry by choosing a gauge. Using the 
isometry (eq. 2) we can first translate the tetrahedron so that one of its vertex, say go, 
is at the identity. This fix one SU (2) invariance, the other corresponds to rotation of the 
tetrahedron around the identity. We can gauge out part of this invariance by fixing the 
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direction of one edge say (10). This still let the freedom to rotate the tetrahedron around 
that edge, we then choose to fix the direction of the plane (120). In term of the variables 
9i,ni the Haar measure reads dgi = ^ sin^ 9id9i(Pni, where (Pn is the normalized measure 
on the 2-sphere. So after fixing go — 1 the measure d/i — dgodgidg2dg3 becomes 

(jjsin'^ 0id0^ d%. (13) 



The residual gauge invariance is fixed by imposing 

/ cos ^13 ^ \ 

sin ^^13 cos 023 
\ sin 6*13 sin 623 / 




n-i^ I U , n2 = sm 6^12 , ^3 




(14) 



In terms of this variables the measure is now 



diJ, = ^ \ rr sin^ 0id0i I sin 0i2d0i2 sin 0i3d0i^ d023- (15) 
\ti J 

It can be checked directly that J d/i = 1. Using relations (9) we can express the measure in 
terms of 0i and 0ij 

2 f ^ \ TT,-^,- sin 0i-id0ii 

This trivially follows from a change of variable and the following fact \{ni A n2) ■ ^3! = 
sin 6*12 sin 6*13 sin 6*23. If we make the change of variable % — >• 0ij, eq.(16) becomes 

^ vr4(nLsm^i)|(niAn2)-n3r 
We can now show that the denominator is indeed 

•>/det[cos 0ij] = ^JJ sin 0^ \ {Hi A 712) • ^3] (18) 

This is done by considering the 4x4 determinant det[cos^/j], which can be written as 
a 3 X 3 determinant det{cos0ij — cos^jcos^^j) by appropriate line substraction. Moreover 
cos 0ij — cos 0i cos 0j = sin 0i sin 0j cos % so det [cos 0ij] = (11^=1 sin^ 0i) det{ni -nj). If we call 
n1 the component of fii in any given basis we see that det(fii-fij) = det(n^)^ = |(ni An2)-n3p. 
Thus we have relation (18). Overall this leads to the measure 

2Y{ jMud0u 
TT^ ^det[cos6/j] 

Domain : The domain of integration arise naturally since, for instance, the relation (8) 
imply that 

cos(^i + 0j) < cos 0ij < cos{0i -0j), (20) 
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which is equivalent to the condition 



\0i - ej\ < Oij < TT - Itt - - Oji- 



Integrand: Finally, we get the theorem using the character formula 

sin(/ + 1)^ 



xiigm 

This completes the proof of the theorem. 



sin^ 



(21) 



(22) 




FIG. 1: Dual angles in a spherical tetrahedron. The dual 623 of 6*23 in the dashed spherical triangle 
appears as the dihedral angle of the edge 01. 



Remarks: 

1) It is interesting to note, by a direct computation, that in the gauge where go = 1 we have 



V'det[cos6'7j] = 1 1 tr( [5(1, 5(2] 5(3) 



(23) 



Also, from relation (18) it is clear that if det[cos ^/j] = then = 0,7ror |(ni Afi2) • ^al = 
in which case the tetrahedra is degenerate. The reverse proposal is trivial. 

2) The integration domain is invariant under the transformation 



f/j CTiajOij H TT, 



(24) 



with a I — ±1. Under this transformation the Graham matrix becomes ajaj cosOu, and its 
determinant is invariant, while the product of sinus produces a factor 



(25) 



KJ 
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The exponent is always a even integer due to the admissibihty condition that the label 
around a vertex must sum into an even integer. Using this extra discrete symmetry we can 
restrict furthermore the integration domain to the domain which is the subset of T>„ such 
that 0i < 7r/2. 

3) In this remark we clarify the fact that we will consider the variables 9ij both as lengths 
and dihedral angles of spherical tetrahedra. There is a duality between points and 2-spheres 
in S^, namely if (? G SU{2) we can define the dual 2-sphere Sg = {x ^ SU(2), tr(a;5f) = 0}, 
Sg is the sphere of radius 71/ 2 centered at g. Reciprocally, such a sphere determines 2 points; 
its center g and its antipodal center —g. So, given a set of angles 9ij G V^,^ we can construct 
-modulo gauge invariance- 4 group elements gj G SU(2) such that cos Ojj — {l/2)tT{gjgj^). 
We can consider this group elements as being the vertices of a spherical tetrahedron, denoted 
T{6), so that 61J is the length of the edge (/J). We can also consider these group elements 
as defining 4 spheres forming the 4 faces of a dual tetrahedron denoted T*{6). In this case 
9ij is the dihedral angle of the edge dual to (IJ) of the tetrahedron T*{9). But given a 

tetrahedron T{9) we can also compute the dihedral angles Ou of this tetrahedron as in 9, 

and construct the tetrahedron T{6ij) : it turns out that we do not get a new tetrahedron 

since T{9jj) = T*{6jj). It can be checked that the Graham matrix of the dual tetrahedron 

T*{9jj) = T{9ij) defined by Gu — cosOu is related as follows to the Graham matrix 
Gij — cos 9ij. Let 9jj be in the interior of then G is invertible, (G"^)// > and 



Gij = , (26) 

^{G-^)n{G-^)jj 

Finally if 9ij is in the boundary of P^, T[9) becomes a degenerate spherical tetrahedron 
and T*{9) is a fiat tetrahedron. 



B. Analysis of the integral 

Our aim is now to study the large N asymptotic of the integral expressing the 6j-symbol 
in terms of invariant variables 

The form of the integrand - an oscillatory function with an argument proportional to N 
times a function independent of N- suggest that the asymptotic behavior can be studied by 
stationary phase methods. However, this is not so simple since the integrand is singular on 
the boundary of the integration domain. We therefore have to make the asymptotic analysis 
by carefully taking into account this singular behavior. For that purpose, wc will split the 
integration domain into two different parts, which will be analyzed by adapted methods. 
We split the integration domain into 

V, = V<^ U with :t _jy _ jyl (28) 

with e > is sufficiently small. Thus the integral splits into and 

respectively given by the integrations on and "D^^. Our main theorem summarizing the 
results obtained for these asymptotic analysis is the following 
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Theorem 2 Let FT*{lij) be the flat tetrahedron which is such that the length of the edge 
(IJ) is given by Ikl (with I,J,K,L all distinct). Let Qrl be the angles between outward 
normals at edge (IJ) , and denote by V*{lij) its volume. Suppose that V*{1) ^ and that 
< e < min(6/j, tt — 6/j). Then the first terms in the asymptotic expansion of both integrals 
are 

T>(N1 \ ^HT.i<j{Nlij + 

^ ^^^''^ ?,^NW*{kj) ' ^^^^ 



and 



Overall leads to the asymptotic for the -symbol 



^(^'"'-ivJl^-^IE^e-^l)- (31) 



This result is, of course, consistent with the results of Ponzano-Regge and Roberts. The 
methods are however very different, and to our taste simpler in the sense that they can be 
extrapolate to the Lorentzian case as we will see in section III. 



C. Asymptotic of />(iV//j) 

Method : The result (29) is proved by the stationary phase method. We recall first the 
main steps of this method on a one dimensional example. First we will write the integral 
under a form which is appropriate for the stationary phase, namely something like 

In= dx (i){x)e'^^^''\ (32) 



with / an analytic function. This is done by proposition 1. The asymptotic expansion of 
such an integral is given by contributions around the stationary points of the phase, i.e. 
points in the complex plane such that /'(ajg) = 0. In our case such points are given by 
proposition 2. The contribution of such a point is then obtained by expanding the phase 
around it at second order and extend the integration to infinity 



/+00 
diSx) 0(xo)e^^(^(^°)+^^"(^°)(*^)') (33) 



0(x) is the part of the integrand independant of N and is evaluated on the stationnary 
points xq. The constant term can be factored out of the integral and the integration is now 
a gaussian which can be performed 



In 



-^5:,(..)e*-r,«„ (34) 

Note that if xq lies on the boundary of the domain, its contribution has to be half counted. 
For the n-dimensional gaussian integration (see eq. 35), we need to compute the determinant 
and the signature of the quadratic form obtained by the expansion. All of this is done in 
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proposition 3. Recall that in the n dimensional case, the gaussian integration is performed 
using the fact that for a real n x n symmetric invertible matrix A with signature cr{A) we 
have 



[dXi] exp 



n" 



detAl ' 



(35) 



Proposition 1 The integral {Nlu) can be written as 



with 



{eij=±} \I<J 

^iT^KjeijOij QiNij^j^jeijlijeij+Y.i^jXicoseijXj) ^-^ 



(36) 



(37) 



Proof : — 

We first split the sinuses into exponentials, which leads to integrals which are still convergent. 
This is shown by notice that due to the computation of the measure, the integral 

[dOu] J „ . (38) 
•>/det[cos6'/jJ 

is just the integral 

dhidh2dhs 



(n.sin^,)(n.<.sin^.,) ^^^^ 

As the domain is a compact one, and exclude the points 9i,9ij = 0,7r, this integral is 
convergent. Now we use the equation (35), to rewrite the denominator as arising from such 
an integration (the signature of the matrix being 4) . One finally rescales the Xj by y/N. 



We consider now each of the integrals and apply the stationary phase method to com- 
pute their dominant asymptotic contribution. The following proposition gives the stationary 

points of the phase for these integrals. They arc expressed in terms of the geometrical ele- 
ments of the tetrahedron FT*{lij). In the following we denote V* its volume, Aj its areas 
and 0/j the angles between its outward normals. 

Proposition 2 The stationarity equations of the integral are given by 

J^cosejjXj = 0, (40) 

2XiXjsm9jj = euhj. (41) 

These equations admits solutions only when ejj = {—l)'^ai(7j with aj = ±1 and r] — 0,1. 
For each of the 16 admissible choices of eu there are 2 solutions s = ±1 given by 

e,j = <J,ajBu + ^"^""■' tt, (42) 

X, = «V,-|i=. (43) 
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Proof : — 

Localisation of the solutions in the complex plane: We begin by looking for general solu- 
tions in the complex plane. Indeed in the method of stationary phase, one should look for 
stationary points in the whole complex plane, as the contours of integration can be analyti- 
cally deformed in C. We prove that a general solution {9ij, Xj) of the stationarity equations 
is in fact such that 9ij e M and e M or iR. We start from the following lemma 

Lemma: Consider ojj a 4 x 4 symmetric matrix such that Ylj — 0- Consider the 5x5 
matrix K[a.) defined by 

KKL^^{aiiOijj-{aijf), Ki.^l, ^55 = 0, (44) 

then we have 

{Coiijf = SiauY detK (45) 

where Cof/j denote the (/ J)-cofactor of K. 

We apply this lemma with ajj = XjXj cos 9ij (the hypothesis is satisfied due to the 
first stationary equation). In that case K is just the Cayley matrix C/j = (due to the 
second stationary equation), whose cofactors A/j are real. As detC > 0, eq.45 says that 
aij = XjXj cos 6ij is real. XjXj sin Ojj is also real from the second equation. Summing 
the square of these two results, 

< {XiXj sin Oij? + {XjXj cos Ou? = {XiXjf (46) 

we can conclude that XiXj is real, and thus that du is real. Lets denote Xi — e^^\Xi\. If 
we look at the phase of the second equation (41) one gets that 

e'^'^'^' = ±eij (47) 

This imply that 

e^^^^= (^^y = L (48) 

Therefore one has e^^' = i^aj with o"/ = ±, 77 = 0, 1 and e/j = (— 1)''o"/(Tj. The conclusions 
of our study of the localization of the solutions are that the only integrals with stationary 
points are the 16 configurations such that e/j = {—l)^aiaj and that the solutions Ojj are 
real, while Xi e i''aiR+. 

Explicit solutions: Now we compute the explicit form of the solutions. Consider a solution 
for a admissible configuration ejj — {—l)^a][aj. The solutions of this configuration and of 
the configuration all e/j — +1 are mapped to each others by 

X'j = iV/Xj (49) 

Thus we only need to study the solution for the configuration e/j = 1, the solutions for the 
others configuration will be recovered by this map. 

The matrix Gjj = cos 6ij can be written as arising from the scalar product of four unit 
vectors nj of M"^, i-e Gjj = nj ■ fij. The first relation (40) imply that det(G) = so the fij 
span a vector space of dimension at most three. It is easy to see that under the hypothesis 
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of the theorem it is exactly of dimension three. Lets denote Xj = XjUj, the second relation 
imply that |X/ A Xj\ = The LHS represent the area of the triangle If the 

vector space was of dimension two or less then these area, hence //j, would be related by 
triangular equalities. This is not allowed because the tetrahedron FT*{lij) is supposed to 
be non-degenerate. Therefore the matrix G as exactly one null eigenvector. 

Thus [11], this matrix G is the Graham matrix of a flat non-degenerate tetrahedron, 
whose angles between outward normals are the 6ij and proportional to the unique 

eigenvector, the Xj in our case. This tetrahedron is defined up to scale. To fix the ideas, 
let us consider such a tetrahedron of volume 1, let us call A/j its lengths, and aj its areas, 
we thus have Xj — iiaj. It is well known (see appendix A) that these geometrical elements 
satisfy 

^^cos9ijaj = Q, (50) 
J 

aiajsmOij = (3/2) A/j. (51) 

By considering this together with the second equation 

lij = 2^jL^aiaj sin 9ij = S/ZXij, (52) 

we see that this tetrahedron is in fact homotethic to the tetrahedron FT*{lij). The Ou 
solutions arc thus ©jj the angles between outward normals of FT*{ljj). As its geometrical 
elements satisfy 

AiAjsmQu^lluV* (53) 

we have 

Xi = s^=Ai, s = ±1 (54) 

Finally for this configuration e/j = +1, there are two solutions. According to the map (49), 
for each admissible configuration this gives the two solutions given in the proposition . 



Our task is now to expand the phase around the solutions at the second order, and 
compute the determinant and the signature of the corresponding quadratic form. 

Proposition 3 For each of the 16 admissible configurations eu — {—l)^aiaj the expansion 
of the phase 

J2 ^uhjOij + 5^ ^/ cos OijXj (55) 

KJ I,J 

around any of the two s = ±1 solutions gives 

KJ 

where Q is a quadratic form. Its signature is a{Q) ~ 2 and its determinant is detQ = (|V^*)^. 
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Proof : — 

Consider first the expansion for e/j = +1, the constant term is easily seen to be Xl/<j hj^u- 
The computation of expansion with fluctuations {59ij, 6Xi) at the second order leaves us 
with a quadratic form, which can be expressed as such 



KJ 



i^j 



E 



K 



1 



Aj cos Q IK ' 



{6Xjf 



cos 0/ J 



5Xi5Xj 



(57) 



where the 6Pij are a shift redefltion of the 69ij 

Consider now the expansion of the phase of an integral le around the corresponding 
solution. One has just to carefully add the aj and rj into the previous results. The constant 
term becomes 



€ijlijiaiaj9ij + 

KJ 



[1 - aiaj) 



(yjfyj 



(58) 



KJ 



KJ 



The second term is always mod. 27r, due to the condition that the labels around a face 
must sum into an even integer. The quadratic form is unchanged, except one has just to 
redefine 8Xi ajSXj to get the same expression. 

Our task is now to compute the determinant and the signature of the quadratic form Q. 
This quadratic form is given by a diagonal part (involving the S(3ij), whose determinant is 



(nj<j^Mjcosgjj) 
(3T/*)6 



(59) 



and a non-diagonal part which is the 4x4 quadratic form involving only the SXj. This 
non-diagonal part can by computed in terms of the Ar and cos 9ij. Using the formula (see 
the appendix A) 



Aj A J cos 9 IJ = -^^ij 
16 

/ 16 ' 



(60) 
(61) 



(A/ J denotes the (/, J) cofactor of the Cayley matrix), the total determinant can be rewritten 

as 



detQ 



1 



(3T/*)6 



Pi^ij) 



(62) 



with P an homogeneous polynomial of order 6. The A/j can now be expressed in terms of 
the lij, the same can be done for the volume V* and symbolic computation (Maple) allows 
to check that 

PiAij) = (63) 
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Overall leads to 



detg ^(Iv*] (64) 



The proposition concerning the signature of the quadratic form is verified only numer- 
ically. We have computed numerically the signature of the quadratic form Q for many 
randomly chosen tetrahedra. In all cases we obtain that a{Q) — 2. 



We can now put together the results of the previous propositions to prove the asymptotics 
of the integrals /g, then of the total integral />(M/j). Consider an integral 7^, if the ejj can 
not be written e/j = {—l)^aiaj, the integral has no stationnary points and its asymptotics 

is a o{N~^^'^). If ejj = (—l)^ajaj, the oscillatory phase has two stationnary points given by 
proposition 2. Recall that the integrand of contains a factor indcpcndant of N (see eq.37) 
which has to be evaluated on the stationnary points. This evaluation is 



exp 



^e/j (ajajQij + (1 - o'/crj)|j 



I 

. KJ 




^e^(-i)''E.<.e.. TT^^ \ (65) 



Recall that these points are the dihedral angles of a flat tetrahedron and are on the boundary 
of the domain according to the remark made after theorem 1. Their contributions have to 
be half-counted. Using the results on the determinant and signature of Q, we obtain that 
the asymptotics of /g is 

^^(^^^^) - ^ (n^^^) ^^e^^-)''S.<.(^^..+i)0.. (66) 

Now there are 8 configurations of the e/j with (—1)'' = +1 and 8 with (—1)'' = —1. Summing 
this according to the expansion (eq. 36), the total contribution is 

which proves the first part (29) of the theorem 2. 



D. Asymptotics of I<{Nlij) 

Method : Before going on with the proof of the second part (30) of the theorem let us 
illustrate the problematic on a one dimensional example. Consider the simple case of a one 
dimensional integral : 



b 

iNx 



dx /(x)e*^^ (68) 



If / is a regular function on [a; 6], this integral has a 0{j^) asymptotic which has its origin 
in the contributions of the points a, h of the lying on the boundary. Now lets consider a 
function which is regular in ]a; h[ but posses a singularity on the boundary, for instance 

\. J^J^', (69) 



a 
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with / regular and non- vanishing at a. one can show that the asymptotic is an 0{1/N^~'^). 
This suggests that the leading asymptotic behavior of an integral with singular points come 
from the contributions around the most singular points. We expect that in general the 
dominant contributions will come from singular points of the integrant. It is easy to show 
that the zeros of highest order of the denominator det[cos^/j], which will dominate the 
large N behavior of the integral I^{Nljj), is given by the configuration where all 6ij are 
zero. The proof amounts to write the integral as a one dimensional integral which realize 
an expansion around this singular configuration. We first express the integral as a one- 
dimensional oscillatory integral over an auxiliary variable S. This is done in proposition 4. 
Then we compute the large N asymptotics of this integral in proposition 5. The result is 
thus expressed as an integral over Euclidian tctrahcdra, involving the inverse of the volume 
V. This integral is exactly computed due to a remarkable self-duality property given by 
theorem 3. 

First it is useful to note that if e is small enough (i-c e < 7r/3) the domain V^^ splits 
into eight disconnected components which can be mapped into each other using the discrete 
symmetry (eq. 24). Using this property wc have 

\7^J J A ^/det[cos0IJ] t<j 
where T>^ is the subset of [0; e]^ of Ojj satisfying the relations : 

0IJ < 0IK + 0JK (71) 
for any triple (/, J, K) of distinct elements. We can now rewrite the integral. 
Proposition 4 

r+oo 

I<{Nlij) = / e'''^S^Fi{S)dS (72) 



where Fi{S) is a function of S depending on the lengths Ijj and given by 

Fi{S) = -^,v.p I E - E n dOu, (73) 

,52./detcos^7j 



eu=±i I J -5V[detcos^7j] 

where v.p means that we take the principal value of the integral i-e v.p J d9 = linia^o J^^^ dO 
and 5{x) is the Dirac functional. 

The rewriting of the integral is immediate, the key point here is that the function Fi[S) 
is well defined since the functional S{9ij) = ^uhjOij don not possess any stationary 
points in the domain "D^ as was shown in the section II C. 

The asymptotic expansion of the resulting one-dimensional integral over S can now be 
computed. 



Proposition 5 



where V{uij) is the volume of the flat tetrahedron which is such that the length of the edge 
{IJ) is given by ujj. The domain is the set of all Euclidian tetrahedron. 
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Proof : — 

It is well known that in the case of a one dimensional integral the asymptotic is: 

/ + 00 f + OO 

dSe'^'^S^FiiS) ~ / dSe'^'^S^Fi{Q). (75) 
-oo J —oo 

Thus what we need is an evaluation at 5" = of the function Fi{S). This function indeed 
regular at this point, as it can be seen by changing Ojj — Sujj. We get 

As yS" — > 0, the domain becomes the domain of all Euclidian tetrahedra. The denominator 
can now be expanded as 5 ^ 0. The proposition follows from this expansion 

^det[cos{Suij)] -5^0 S^6V(uij) (77) 

where V{uij) denotes the volume of tetrahedron, whose length (/J) is ujj. This lemma 
can be proved by direct expansion, but if follows from this observation : Consider the 
determinant A = detcos^/j (with ^oi — ^i)- By substracting row 4 (resp. hne 4) to the 
three others, one obtains the determinant 

A = det (cos 9ij — cos 9i cos 9j ) (78) 

which is rewritten, following relation (eq. 11), 

A = ( sin^ ei)det{ni ■ Uj) (79) 

i 

The limit — > of this expression leads to the expression det{9i9j{ni ■ Hj)) which is the 
square of the volume of the flat tetrahedron defined by vectors 9ini, ^2^2, ^3^3 (see eq.A9). 
Thus we have the property 

det[cos^/j] ~e^o eV^ (80) 
with V the volume of the tetrahedron whose lengths are {9i, 9ij). 



We have express the asymptotic of the integral I'^{Nlij) as an integral (74). This integral 
can be in fact exactly computed, showing a self-duality property of the inverse of the volume 
of the tetrahedron. 

Theorem 3 Let V{uij) (resp. V*{lij)) he the volume of the flat tetrahedron which is such 
that the length of the edge {I. J) is given by uu (resp. Ikl) The inverse volume satisfy the 
remarkable self duality property 

A' f , HkjMIijUij) 1 

The proof of this theorem make the use of the following proposition concerning the 
measures on the spaces of Euclidian tetrahedra. 
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Proposition 6 Let Ui be 3 vectors ofW^, and lets denote uoi — \ui\, Uij — \ui — Uj\ then 

. 27r ?>'kV{uij) ^ ' 

where V{uij) is the volume of the tetrahedron {ui,U2,U3) whose edges length are ujj. 



Proof : — 

The measure on is given by d^Ui — ATrufduid^fii, where Ui — \ui\, fii & and d^n is the 
normahzed measure on the 2-sphere. So the measure d/j, — Yli becomes 

dfx = 2^ (jj u'^idu^ d^fii. (83) 

Prom there we proceed exactly as in the proof of theorem 1, using the fact that 

^V{uij) = (Y[ui)det{ni A ■ ns). (84) 



We can now prove the theorem. This proposition imphes that the integral eq. 74 has a form 
similar to eq. 1 

2\^ 1 



where G is replaced by and the character Xiio) is replaced by Ki{u) = sin(/|-u|)/|-u|. This 
fact was already conjectured in the paper of Baez et al [8] where the kernel Ki{u) has been 
interpreted in term of spin networks of the Euclidean group. Changing the lengths to their 
dual li — Ljk, lij — Lfe and using the Kirillov formula 



sin L 


\x\ 




\x\ 





e'^^-"ci"n (86) 

52 



one obtains 



I'^ihj) = AdT ^ij) I TT d^'^i I (TT d'^nij) TT exp {iuk ■ {UjUij + LjUj - Ufii)} 
i<j -^^^ hs^ff^j Y 

(87) 

The integrals over the sphere can be rewritten as integrals over R^: 



/ 



I d^X5{\X\-L), (88) 



while the integrals over the Ui can be performed J e^'^^d^u = {2ti)^6^{x), leading to others 
5-f unctions. 

(89) 



15 



The integrations on X^j can be performed, due to the 5^^^ functions. It remains 

I^ilij) I d^iU - Li) n T^^d^^ - ^il - Lij) (90) 

— * — * — * 

Changing into variables Xi — \Xi\,Xij — \Xi — Xj\ and using again the proposition 6 one 
gets finally 

'"^"""^ ~ 3.]V»V(L,L,) = SWV5F(W 

Remark : The self duality property (81) of the inverse volume shares strong similarity 
with the self duality property of the square of the quantum 6j symbol which has been 
shown by Barrett [12]. It is therefore tempting to ask whether this property has a natural 
interpretation in terms of quantum groups where the square root of the inverse volume would 
be interpreted as a quantum 6j symbol. 



III. LORENTZIAN 6J SYMBOL 

In [2] it was shown that the partition function and transition amplitudes of Lorentzian 
3d gravity can be written as a state sum model. To construct this state sum model one 
triangulates a 3d manifold, one gives an orientation to the edges of the triangulation and 
color these edges by unitary representation of SL(2,M). The weight associated with an 
oriented colored tetrahedron is the corresponding 6j symbol of SL(2,]R). Among all the 
possible representations of SL(2, IR) only the positive discrete series, the negative discrete 
series and the principal series play a role. All these representations are unitary and infinite 
dimensional [13]. We will denote Tji+, / G N, the positive discrete series of weight /, T_ii-, 
I e N, the negative discrete series of weight /, and Tp, p G M"*", the principal series of weight 
p. To simplify the exposition we have adopted a unifying notation, all the representations 
are denoted Tx where A = or p refers respectively to positive discrete series, neg- 

ative discrete series and principal series. Similarly to the case of SU(2) one can construct 
the Clebsch-Gordan coefficients, and 6j-symbols as a recombination of Clebsch-Gordan co- 
efficients (we refer the reader to [14] for a precise definition of these objects). Due to the 
fact that we can color each oriented edge of the tetrahedron by either il'^ , or p there 
are many different types of 6j symbols. However the square of all these 6j symbols can be 
expressed by an integral formula similar to eq.l 



^(A/j) 



I Aoi Ao2 Ao3 I ^ ^ / Ux>^.A9J97')[Ud9 

{ A23 Ai3 A12 J J(SL(2,M))4 f_^j V 



(92) 



where X\{9) is the character corresponding to the A representation (see appendix B). In 
the following we will suppose that the tetrahedra is ordered by the vertices order, i-e (/J) 
is positively oriented if / < J. The triples A/j, Xjk, ^il are supposed to satisfy some 
admissibility conditions described in the appendix B. The measure [Y[j dgi] denotes the, 
fully gauge fixed, product of Haar measures. It is here necessary to fully gauge fix the 
symmetry gi — > kgjh in order to avoid divergences since the group is non compact. Using 
techniques similar to the one we used in the previous section we can compute the gauge 
fixed measure 
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Proposition 7 



[TT dgi] = ( n dm] (93) 

^ ^ / V|det[G,,(T)]| ^ ' 



where tjj — iOu and Gjj — cos 9jj if gjgj^ is elliptic and conjugated to ke^j in this case 
drij = dOij; and tjj = G/j = uij coshtu if gjgj^ is hyperbolic and conjugated to 

vijatjj, fij — ±1, see appendix B for the notations. 

These gauge fixed variables are interpreted as parametrizing the invariant geometry of an 
AdS tetrahedron. The technics used to deal with non compact spin network integrals where 
first introduced in [15]. Similarly to the case of SU(2) one can split the integral expression 
of the Lorentzian 6j into a sum /^(A/j) + /<(A/j) where /^restrict the integration range 
into a domain where |r/j| > e and \iT: — tjj\ > e. In the following we analyze the asymptotic 
behavior of I^{N\jj) when oo. In order to illustrate our method, we will treat in 

detail the case where all A/j are real, corresponding to a spacelike Lorentzian tetrahedron, 
and the case where all A/j are imaginary, corresponding to timelike Lorentzian tetrahedron. 
The last subsection gives a method to solve the stationnarity equations in the general case. 



A. Spacelike tetrahedra 

The case where all representations are continuous is simpler since, first there is no restric- 
tion on the set pu and second the characters Xp ^-^e supported only on hyperbolic elements. 
The integral reads 

i'ipu)-y:f -^^^^^i^^iUdtu]. (94) 

„jj=±Jtij>e V|det[i//jcosht7j]| 

We have uu — +1 and ujj are restricted to the ones actually arising from an AdS tetrahe- 
dron. The cosine can be expanded and the denominator expressed as a Gaussian integral, 
overall this leads to 

where the oscillating phase is 

^XiVijCOshtijXj. (96) 



KJ I,J 



The asymptotic behavior of this integral is driven by the stationary points of this action. 
They satisfy 

J2^ijcoshtijXj = 0, (97) 



J 



2XiXjsmhtij = -^ijeupij. (98) 

The analysis of this system of equation is similar to the one performed in proposition 2. This 
analysis will be done in subsection III C in the general case of a Lorentzian tetrahedron. The 
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conclusions are given in terms of the geometrical elements of the tetrahedron T(p) given 
by the spacelike lengths pij. Let us denote Ai its areas, V* its volume and {Vij,Tij) its 
Lorentzian dihedral angles. The possible solutions are given by 

X, = ±^V.-|= (99) 
tij - Tjj (100) 
The expansion of the phase around the stationnary point T/j gives the constant term 

-5 = ^e,jp7jT,j (101) 

KJ 

The quadratic form obtained by expansion around these configurations has a determinant 
which equals to — (|V^*)^- In order to prove the result about the determinant one first 
write the quadratic form in term of cofactors A/j of the Cayley matrix (see appendix A). 
The form we obtain is the same as the one obtained in the SU(2) case except for a global 
minus sign in front of the non-diagonal part. It was shown there that the determinant 
of the quadratic form is proportional to the fourth power of determinant of the Cayley 
matrix. Since this is a polynomial identity it extend when the Cayley matrix comes from a 
Lorentzian tetrahedron. Numerical simulations suggest that the signature of this quadratic 
form is ±2 depending on the number of positive and negative Vjj in the dihedral angles of 
the spacelike tetrahedron T(p). We now have to consider the compatibility equations on 
the signs to actually determine the possible cri,i'ij, eij,r). The first stationnarity equation 
can be written 

Ajajuij cosh Tij = ^ AjajUijVijGu = (102) 



E 



where Gjj = Vjj cosh Tjj is the Graham matrix of the tetrahedron T{p) , while the sign of 
the second equation gives 

eij = -(-l)V,j(7j(7j (103) 

Considering the equation (102), we know (see appendix A) that the Vjj can be written 
Vjj — —ajaj. If we consider the factorizable configurations of i/jj — PiPj, we get for the 
equation (102) 

-aiPi Ajajaj(5jGij = (104) 
J 

The fact that Gu possess only one null eigenvector given by the Aj amounts to identify 
Pj = ajaj and the contributing uu as uu = —ajajVij. Considering now the second 
compatibility equation (103) gives the possible e/j as ejj — Taking into account 

theses contributions, we obtain the following asymptotic term 

, C MEkj ViJPijTij) 
for C a constant independant of the pij. 



18 



B. Timelike tetrahedra 



In the case where the representations are discrete and equal to ^ the integral reads 



I>{iNeijr/f) = / n ' ' . ' ^[TT dOij] + ■■■ (106) 

where the domain of integration is the set of angles 9jj G [e, tt — e] U [— vr + e, — e] which 
are such that det [cos 6'7j] is negative. It corresponds to the set of timelike AdS tetrahedra. 
The dots refer to the fact there are other terms in the integral coming from the sectors were 
at least one element gjgj^ is hyperbolic. These terms contain, in the integrand, at least 
one factor exp— A^|t/j|, such a factor is bounded by cxp — A^e and is exponentially small, 
therefore they do not contribute to the asymptotic behavior of the integral. Even if the 
expression eq.106 is suitable for the asymptotic analysis, it is interesting to note that we 
could restrict the integration to be the set of angles for which 6ij e [e, tt — e], det[cos 9ij] < 
if we use the identity 

Pie y \ de^ / p(e)sm(Ke) de. (io7) 

smB 2i Jq 

for P such that P{—0) = P{0), which is the case of the denominator for each of its arguments 
9jj. The integral now reads 

a form which is strikingly similar to the Euclidean integral eq.4. 

The analysis of the asymptotic behavior goes along the same line as before and the 
stationary points are in one to one correspondence with timelike tetrahedra. Modulo per- 
mutation of the vertices and change of orientation one can always present such a tetrahedron 
in a form where all the edges are future timelike vectors and the vertex 3 is in the future of 
2 which is in the future of 2... in the future of 0. This amount to take al A/j, / > J to label 
positive discrete representation. For such a configuration and using again the same technics 
one expect the asymptotic to be 

n^NliJ) ^ + C.C (109) 

where 0/j are the dihedral angles of the timelike tetrahedron whose edge are //j, V* is its 
volume. C is a constant independent of the Ijj and c.c stands for complex conjugate. 



C. General method for the stationarity equations 

In the general case of Lorentzian 6j symbol, the edges are labelled by (discrete or con- 
tinuous) representations. If the representation is continuous then the character is zero on 
elliptic elements and if the representation is discrete the character is exponentially small 
on hyperbolic elements. Thus the only case of interest for the asymptotic behavior is when 



19 



the integration is over elliptic elements if the edge is labelled by a discrete representation 
and over hyperbolic elements if the edge is labelled by a continuous representation. This 
amounts to integrate over AdS tetrahedra which have spacelike edge when the representa- 
tion is continuous and timelike edge when the representation is discrete. The asymptotic 
behavior of this integral is controlled by the solutions of the following stationary equations. 

Y,XjVijcoshTij = (110) 
J 

2XiXjUijsmhTij = -eijXij, (111) 

In the elliptic case, A/j stands for ±ilij, tjj denotes i9ij and ujj is always 1. In the 
hyperbohc case Xjj = pu, tjj = tu and ujj — ±1 for the different sectors of hyperbolic 
angles. We will express the solutions of these equations in terms of the geometrical elements 
of the Lorentzian tetrahedra T(A) given by the square lengths Xjj. Let us denote by 0/j and 
{Vij, Tij) its dihedral angles respectively in the timelike and spacelike case. Let us denote 
Aj its areas and V its volume. We will procede as in the SU(2) case by first identifying 
where the solutions lie in the complex plane, then compute them explicitely. 
Localization: We can repeat the analysis done in the SU(2) case. We apply the same lemma 
(eq.45) to obtain 

with C(A) the Cayley matrix of the Xjj and Ajj its cofactor. This cofactor is real and the 
determinant of C is negative. We can conclude that XjXj cosh tjj is real. Now the second 
stationnarity equation tells us that XiXjsmhru is real in the spacelike case and pure 
imaginary in the timelike case. Taking the difference of the square of these two equations 
gives 

{XjXj)' = {XjXj cosh Tij)' - {XiXj sinh nj)' = _g'^^^c{X) ~ ^ ^^^^^ 

One can conclude that XjXj is real since the RHS of this equation is always positive. This 
is clear in the timelike case since the square of Xjj = iljj is negative. In the spacelike case, 
this arise from rewriting it using the relations (A5) and (A8) 

(A„)= Xh _ ^MJ ^„ _ m r„ ^ ^ (114) 



-8detC(A) 4 (3^)2 (3^)2 {3Vy 

We thus know that XjXj is real, cosh tjj is real, sinh r/j is real in the spacelike case and 
pure imaginary in the timelike case. This allows to conclude that t/j is pure imaginary in 
the timelike case (hence 9ij is real) and t/j is real modulo in in the spacelike case. However 
this in ambiguity corresponds to a ± ambiguity in the sign of cosh r/j and sinh r/j, and 
thus leads to the same solution obtained by taking tu and changing the sign of uij. These 
solutions are already taken into account and we can restrict to the case of r/j real in the 
spacehke case. This completes our investigation of the locahzation of the solutions. 
Explicit solutions: We consider first the equation (112) ior I — J and use the equations (A8) 
and (A7) to rewrite the RHS. We obtain 
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Thus it exists aj = ±1 and rjj = 0,1 such that Xj = aii'^^ /y/SV. The fact that XjXj is 
real for all IJ leads to the fact that all rji are equal and 

Xi = aii^^L (116) 



Now we consider the equation (112) for 7 J and use the solution we found for the Xj 
to get 

(coseuY = {cosOuT (117) 
{coshtijf = (coshrrj)^ (118) 

We have to solve these equations in R, but since the original integral is symmetric by 
changing t — > — t and 9^9 + 27r, we keep only the actually different solutions and we solve 
these equations for t e R"^ and 9 e [— tt; tt]. This gives the solutions 

tij = Tij (119) 

9ij = aij [(5ijQij + (1 - /?/j)|) for au, (3ij = ±1 (120) 

Reporting these result in the first stationnaity equation gives 

We still have to determine the possible uu, uj, Pij... by considering compatibility equations. 
In the previous one, one can use the fact that by hypothesis, the Graham matrix of T(A) 

«'-{rx'"} (^^^' 

possess only one null eigenvector given by the Aj. On the other hand, the sign of the second 
stationnary equation gives also equations 

which should be satisfy in order for the solution to exists. The bottom line is that the 
system of equation 110 admits solutions only if A/j could be interpreted as a set of length 
of an oriented tetrahedra T{\jj). In that case the on-shell action for this solutions is, up to 
a sign, proportional to the Reggc action of Lorentzian tetrahedra. The determinant of the 
quadratic form obtained by looking at fluctuation around these configuration is proportional 
to the volume of the tetrahedron. 



IV. lOJ SYMBOL 

A. Integral expression for the lOj-symbol 

To analyze the case of the lOj-symbol along the same lines of the 6j-symbol, we will use 
higher- dimensional generalizations of propositions used in the previous parts. The analog 
of the theorem expressing the symbol as an integral over invariant variables is 
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Theorem 4 The lOj-symbol can be expressed as the following integral 



^ / ( n + 1)^^-^) ) '^(G'(M) [H (124) 



where G{9) denotes the determinant of the 5x5 Graham matrix Gjj = cosOjj associated 
with the spherical 4-simplex 01234 and S is the Dirac functional. V'^ is the set of all spherical 
4-simplices. 



Proof : — 

We start from the Barrett's expression for the lOj-symbol [9] 

/ XkA9i9l')[Ud9i] (125) 

Then we apply the same method as in the 6j case. By gauge fixing go — 1 we are left with 

an integral over SU(2)^ with an ArfG- invar iance. We have seen in the proof of theorem 1 
that if gi,g2 are two groups elements we can express their invariant measure in term of the 
angles 6^1, 6^2, Ou. 

2 

dgidg2 = sin OidOi sin 62d62 sin 6i2d6i2 (126) 
Moreover if gi is an additional group element its measure is given by 

1 sin 9id9i sin 9ijd9iism92id92i , 

= ^ 

where G\2i is the determinant of the Graham matrix for the spherical tetrahedron whose 
vertices arc l,gi,g2,gi. using these results for g3,g4 allows to prove the following integral 
expression for the lOj-symbol 



lljsin{lij+l)9ij sinjlu + 1)^34 ,^28) 
^"^ VA33 A44 sin 9u 



where the set T is the set 1 < / < J < 4, (/, J) 7^ (3,4). A/j denotes the cofactor (i-e 
determinant of the minor) of the 5x5 Graham matrix. The equation (128) follows from 
eq. (127) since the minor A33 (resp. A44) is the Graham determinant of the 3-simplex 0124 
(resp. 0123). 

One of the main difference with the case of the 6j is the fact that the angles 9ij are not all 
independent. It is clear from our parameterization that 6*34 is not needed : it is a function of 
the 9 others angles. Indeed, from the original integral expression, the 5 integration elements 
gi e can be considered as 5 unit vectors in M^. They are not independent and the angles 
9jj can therefore be considered as the dihedral angle of an Euclidean 4-simplex. This means 
that the 5 dimensional Graham matrix Gij{9) is degenerate, G{9) = det[G/j(^)] = 0. 

In order to get a symmetric expression of the integrant, let T(9ij) be a spherical 4- 
simplex whose dihedral angle are 9jj, denote by G{9) its Graham determinant, and by A/j 
its cofactor. It is clear that 

^^^^^ ^-2Kjj{9)sin9ij, (129) 



39 



ij 
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where the factor 2 is due to the symmetry of the Graham matrix. When the 4-simplex is 
flat (i-e G=0) and non degenerate, the Graham matrix is of corank one, which means that 
the cofactor matrix is of rank 1 

Ajj = AiiAjj. (130) 

This means that the denominator in the integrand of (128) is |934G'(^^)|/2 and the measure 
is simply 

( n deu)—^-^r^^ = {ll^eu)25{G{9)). (131) 

Sm6'34VA33A44 

KJ, 
(7,J)7^(3,4) 

Using this measure in (128) proves the theorem. 



B. Asymptotics for the lOj symbol 

Following the lines of the method employed for the 6j symbol's study, we split the integral 

/ into a sum / = + where in the angles Ojj are restricted to be in [e, vr — e], 
e > small enough. The asymptotic of the integral I^{Nlij) is governed by stationary 
phase. Given Ijj we denote by T*{lij) the Euclidean 4-simplex such that the area of the 
face opposite to the edge (/J) is given by Z/j, and denote the dihedral angles of this 
tetrahedron. If lu are such that T*{ljj) is not degenerate and e < min(©(Z), tt — Q{1)) then 

T>(N1 \ ^ cos(iV(/,j + l)eij) 

where V{1) has to be determined. This property was proven by Barret and Williams [6]. 
Since the argument of the cosinus is the Regge action, it was a strong evidence that spin 
foam models constructed from the lOj give amplitude of 4d gravity. Note that if T*{1) 
is degenerate, then for e > 0, there is no stationary point in the action and I^{Nlij) — 
o(l/7V9/2). 

However, we have to take into account the asymptotic behavior of First, as for the 
6j, we find that the integration domain of I'^ is split into 16 disconnected components. All 
this components are equal and related by the symmetry transformation 6jj ajcrjOjj + 
(1 — (T/(Tj)7r/2. With this symmetry, we can restrict the angles to be in [0,e]. According to 
our analysis of the 6j the asymptotic behavior of this integral is governed by the singular 
contributions of the integrand, i-e points where dijG{9) — 0. The most singular points are 
indeed the ones where 9ij — 0. The behavior of the Graham determinant around this point 
is given by 

G{Su) r^s^o ^'[4!\/4.(m)]2 (133) 

where V4s('u) is the volume of the Euclidean 4-simplex whose edge lengths are M/j. We 
can perform an analysis which is almost identical to the one of section II D, therefore we 
do not repeat it. This analysis shows that the asymptotic of is given by the following 
contribution 

I'-iNlu) - / \{MliJUijmms{u)f) [\{dujj] (134) 

■^^"^ KJ KJ 
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where the integration is on Voo, the set of euchdian flat 4-simphces. 

By an analysis similar to the one done in the previous section, we can express the measure 
in term of a fully gauge fixed measure 



IJ 1=0 



(135) 



where ujj = \uj — ui\. The asymptotic contribution of /"^ is therefore 
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I ttI^Itt' 



s\n{ lij\uj - ui\ 
\uj - uA 



(136) 



which is the form proposed by Baez et al 



V. DISCUSSION 

As we mention in the introduction, given a triangulation of a 4-dimensional manifold, we 
can construct state sum models (usually called Barrett-Crane models) where the amplitude 
associated with the 4-simplex is given by the lOj symbols. The importance of these models 
lies in the hope that they would give a way to compute transition amplitudes in 4d gravity 
[4]. Also, it was found that the integral form of the lOj [9] lead to the interpretation of these 
weights as Fcynman diagrams [16]. This leads into the striking result that the sum over all 
triangulation of Barret-crane amplitude could be interpreted as a sum over Feynman graph 
of a non local field theory [17]. 

However this model suffers from two diseases. First, the amphtudes calculated in terms 
of lOj symbols are always positive [7] and second, as we have just seen, their asymptotics 
is not dominated by the semi-classical Regge action but it is dominated by degenerate 
configurations. Does this mean that these model are therefore not good candidates for 
quantum gravity amplitudes? In order to answer this, we have to recall why the Barret- 
Crane models are suspected as good candidates for GR amplitudes. The original derivation 
of the model in [4] did not make any reference to the Einstein-Hilbcrt action, however there 
has been work since then linking the Barrett-Crane model with a discretization of the path 
integral for Euclidean 4d gravity [18, 19]. In these works it was argued that the Barret-Crane 
model corresponds to a discretization of a constraint BF model introduced by Plebanski 
[20] , and it was shown that the partition function of this model was related to the gravity 
partition function if one excludes the degenerates metrics from the path integral. However 
in [19] it was argued that the degenerate sector is likely to dominate to path integral since 
the number of degenerate configurations is likely to be more important that the number 
of non degenerate one. A similar fact is well known in statistical mechanics under the 
name 'order by disorder' [21]. It is best ilhistratcd in the XY antiferromagnetic model on a 
Kagome lattice. This model is frustrated and posses a huge degeneracy (the vacua states 
have a macroscopic entropy). Since all the states have the same energy, only the entropy 
of the fluctuation around these states could make a difference between them. When the 
temperature is small, the states which are selected are the soft modes states, (i-e the one 
around which fluctuation in energy is quartic), whereas the hard modes (the ones for which 
the fluctuation in energy is quadratic) are suppressed. In Euclidean gravity we are working 
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with the Einstein-Hilbert action 

S = je'Ae'A F''\A)eijKL. (137) 

where is the frame field / are S0(4) indices and F^^{A) is the curvature of the S0(4) 
spin connection. If we evaluate the action on solutions of the equation of motion we find 
that S — So the on-shell configurations are not distinguished by their action. By a 
similar argument to the one used in statistical mechanics one expect that the dominant 
configuration of the path integral are given by the soft modes in the limit where — > 0. It 
is clear that the fiuctuations of the action around a configuration where the metric (hence 
e) is non degenerate are purely quadratic. The softest modes of S is given by e = 0, and 
one expect that fluctuations around this configuration dominates the path integral in the 
classical limit if one include degenerate metrics. Note that in 3d the fiuctuation around e = 
are quadratic, so in this case, we do not expect the degenerate contributions to dominate 
the non degenerate ones. 

This argument suggests that the disease we have just find in the asymptotic behavior of 
the lOj symbol is in fact a disease of the gravity partition function if one includes degenerate 
metrics. The cure for this disease is then clear, one has to avoid degenerates contributions. 
How this can be implemented in the state sum model is also clear. We have seen that given 
e > the integral defining the lOj symbol could be split int I'^ + where account 
for the dominant degenerate contributions, whereas has a nice semi-classical oscillatory 
behavior. One can therefore propose a model where, instead of taking the lOj symbol as an 
amplitude for the 4-simplex, one takes the truncated lOj symbol This cures at once the 
two diseases mentioned earlier. 

One of the key feature of the Barret-Crane weight was the fact that it could be inter- 
pretated as a Feynman graph of a group field theory. The modification of the weight we 
propose still preserve this key property since can be written as a Feynman integral with 
the propagator 

K,.i{9u92) = Xi{9i92')Ce{9i92^) (138) 
where xi denotes the character of the representation of weight / and Ce(e^) is a cut-off 
function which is if |X| < e and 1 if |X| > e. This let open the possibility of writing 
the sum over triangulations of state sum model constructed with a truncated lOj symbol as 
a Feynman graph expansion. The drawback of the truncated model is the presence of the 
cut-off e which still deserve a physical interpretation. We hope to come back to these issues 
in the near future. 
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APPENDIX A: GEOMETRY OF THE EUCLIDEAN AND LORENTZIAN FLAT 
TETRAHEDRON 

Let e/ G M^, / = 0, 1, 2, 3 be the vertices of a non degenerate tetrahedron, we denote by 
Ai the vectors normal to the face opposite to the vertex /. They are defined as 

^7 = ^(ej CftT eft: cl cl ej) (Al) 
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where /, J, L is an even permutation of 0, 1, 2, 3, and 

(7xA^v)" = r/»%fee«V, (A2) 

where Sabc is the totally antisymmetric tensor with 6123 = 1 and rjab = f]"'^ is the metric + + + 
in the Euclidean case and — h + in the Lorentzian case. In all the following we denote the 
wedge product simply A and the scalar product with respect io 'qhy u - v without refereing 
exphcitly to the dependence of the metric. The vectors Ai clearly satisfy 



0. 



(A3) 



The edge vectors are given by Ijj 
relation 



ej — ei, they satisfy with the area vectors a duality 
A • 4- = W5ij, (A4) 



where V = (/qi A /02) ■ ^03/6 is the oriented volume of the tetrahedron. Moreover if one uses 
the identity {u Av) Aw — r)[{u ■ w)v — {v ■ w)u], where 77 is the signature of the metric, one 
gets 

AiAAj^ V^Ikl (A5) 

where /, J, L is an even permutation of 0, 1, 2, 3. 

The Cayley matrix associated with the tetrahedron is defined as 

/ 12 1 \ 



c 



'■12 



n2 ns 

^23 ^2 1 

^13 ^23 ^3 1 

q q /2 1 



V 



''2 
1 



1 1 



(A6) 



where I — l-l and Zj = Zqi- This matrix encodes a lot about the geometry of the tetrahedron 
as the following identities show 



det(C) = ri2%6Vf, 

Aj-Aj ^ -^Ajj, 
16 

where Ajj is the cofactor matrix of the Cayley matrix. 

The first relation can be proven as follows, we have the identity 



(A7) 
(A8) 



det(/i ■ /,) = det{ltvabl^,) = det{ltfv = v{QV) 



Using the relation. 



li ' li 



The square of the volume is rewritten as 



2^6V 



h2 



2/2 



q + q 



'^12 + ^1 + ^2 



2/2 



-P +p + p -p 4-]^ 4-] 
'is ^ 'l ^ 'S '■23 ' ''2 ' 



+ ^1 + ^3 

2 

23 



_72 >p,p 
•"JS •'2 ^ ''3 



2/2 



(A9) 
(AlO) 

(All) 
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This 3x3 determinant is easily seen to be the Cay ley determinant. This follows from 
substracting line 4 to lines 1,2 and 3, then row 4 to rows 1,2 and 3 in the Cayley determinant 
to finally obtains 

—2lf ~ If ~ I2 ^13 ~ ^1 ~ ^3 H 



^12 ^1 ^2 



23 ''2 

2 ;2 72 72 72 72 — 



72 _ /2 _ 72 72 _ 72 _ 72 

'13 '1 '3 '23 '2 
/2 72 

n '2 




^9*^ ^9 ^'i ^2 ^ 

q 



'3 




1 

1 



(A12) 



which reduces to the 3x3 determinant. 

The relation (A8) is proven by direct computation and the use of the identity 

{u Av) ■ {w A x) = r)[{u ■ w){v ■ x) — {u ■ x){v ■ w)] (A13) 

Lets denote Aj — \J \Aj ■ Ai\ the area of the face /, lets introduce the normal unit vector 

Hj, Ai — Ajfii and lets denote Gu — fij -nj the Graham matrix. The relation (A3) implies 
that 

J2AjGu = 0. (A14) 
J 

Prom relation (A8), we get that for 7, J = 1.. 4, the cof actor A/j are related to the elements 
of the Graham matrix 

Aij = -r)16AiAjGu (A15) 

A consequence of eq.Al4 is that 



1=1 



ij 







If the tetrahedron is Euclidean Gu — cos 6/j, and the relation (A5) imply 

3 

AiAjsinOij = - \V\Ikl- 



(A16) 



(A17) 



In this case the Cayley determinant is positive. 

If the tetrahedron is Lorentzian and time-like (i-e all its edge are time-like Ijj ■ Ijj = —Ijj) 
then the normal vectors are all space- like and Gjj — cos ©7 j and 

AiAjsmeij = ^\V\{±)lKL. (A18) 

where the sign in the RHS is -|- (resp. — ) if the edge Ikl is future (resp. past) pointing. In 
this case the Cayley determinant is negative. Conversely given a matrix Gu = cos 0/j which 
is degenerate (of corank 1) one can associate uniquely to it a pair of Euclidean tetrahedra 
or a unique Lorentzian tetrahedron. 

If the tetrahedron is Lorentzian and space- like (i-e all its edge are space-like lij-ljj = pjj) 
then the normal vectors are all timehke and Gjj — uijCOshTu, T > 0, ujj — ±1 and 
i/ii = — 1. vij is the sign of the scalar product (rij • nj), depending on whether ni and nj 
lie in the same part of the light-cone. The z//j can thus always be written as z//j = —ajOj 
with aj = ±1. We call the pair {uu^Tij) where T > vu = ±1 the dihedral angle of the 
edge (IJ). The relation (A5) imphes 

AiAjsinhTij = l\V\pKL. (A19) 
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APPENDIX B: SOME FACTS ABOUT SL(2,M) 



We denote by Tj;+, / G N the positive discrete series of weight /, T^u-, I G N the negative 
discrete series of weight /, and Tp p & IR"*" the principal series of weight p. In this paper we 
consider only the representations of SL(2,R) which extend to representations of S0(2, 1), 
this essentially means that / G 2N. It is well known that the continuous representation 
are conjugated to itself, whereas the positive discrete series is conjugated to the negative 
discrete series. So if one change the orientation of one edge this amounts to change the label 
A into its complex conjugate A. Let (Ai, A2, A3) be a triple of representations meeting at 
a vertex of the oriented tetrahedron, modulo a change of orientation wc can suppose that 
all the edges are incoming at the vertex. Such a triplet is always admissible if there is at 
least one continuous representation. If there is no continuous representation the admissible 
triplets (for the incoming orientation) are 

(/^, l^, l^) with h>h+ k and h + k + h G 2Z; 

{l^, I2, l^) with ^3 > /i + I2 and h + k + k & 2Z. 

Let g be an element of SL(2, R). We say that g is elliptic if it is conjugated to the matrix 

, f cos ^ sin ^ \ . . 



sin 9 cos 9 

g is said to be hyperbolic if it is conjugated to the matrix 

±«< = ±(oe-')- 

Note that and k^Q are not conjugated to each other in SL(2,M), whereas at and a_t are. 
The characters of the positive discrete series are given by [13] 

1 f,i{i-^)e 



if g is conjugated to ke and 



if g is conjugated to ioj t > 0. The characters of the negative discrete series are conjugate 
to the positive discrete series characters 

Xii+ig) = x-ii-ig) (B5) 

The characters of the principal series are given by 

Xpig) = (B6) 

if g is conjugated to ke and 

, , cos pt 

if g is conjugated to ±at t > 0. 

Since the representations are infinite dimensional the characters of these representations 
are not function but distributions. The previous expressions should be consider as a particu- 
lar representative of the character distribution in terms of a locally integrable function. One 
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has however the freedom to change the representative function as long as they differ from 
the previous one on a set of measure 0. It turns out that we effectively need to take an other 
representative function in order to make sense of the integral (92), since the integral admits 
singularities around gj = gj. We propose the following modification. Lets consider a > 
and lets define Ca{g) to be equal to l/(e* — e~*) if g is conjugated to ±at, < t < a and to 
be equal to otherwise. For the continuous series p lets consider the following distribution 

< x'plf >= l^ma^o [ ixpig) - Ca{g))f{g). (B8) 
Jg 

If / is regular at the identity then it is clear that < x'plf >=< Xp\f > so x'p define the same 
distribution as Xp- Analogously we define X±i± by substracting the divergence at 0. And we 
use these characters in the definition of the integral (92). It is easy to see that with such a 
representative of the character the integral is convergent. 
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